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Abstract 

The process e + e~ — > v e e~W + — > v e e~ud is considered as an example of 
the problems associated with maintaining gauge invariance in matrix elements 
involving unstable particles. It is shown how to construct a matrix element 
that correctly treats width effects for the intermediate unstable W boson and 
that is both SU{2)l and J7(l) e . m . gauge-invariant. SU(2)l gauge-invariance is 
maintained by Laurent expansion in kinematic invariants and U (l) e .m. gauge- 
invariance is enforced by means of a projection operator under which the exact 
matrix element is invariant. 



1 Introduction 



In recent years the difficulty and importance of producing exactly gauge-invariant 
amplitudes for processes involving unstable particles has come to be appreciated. 
Prior to 1991 all calculations of electroweak physics at the Z° resonance, and even 
the definition of Mz, were gauge-dependent. A variety of techniques have been applied 
to try to remove explicit gauge- dependence. Some of these involve the generation of 
gauge-invariant Green's functions, such as self-energies and vertex corrections [[L], 0, [|. 
Such techniques however replace dependence on, say, the gauge parameter by some 
arbitrary choice in the definition of the Green's function. In the case of the so-called 
'pinch technique' the procedure used may be rigorously justified by the background 
field method @, |j for a particular choice of the gauge parameter, £. The fact that a 
particular choice of gauge is involved indicates that the gauge dependence has been 
disguised but not eliminated. The appearance of gauge-dependent amplitudes is a 
symptom of a perturbative expansion that does not reflect the physics of the problem. 

The only method to produce genuinely gauge-invariant amplitudes at all orders 
in perturbation theory is that of Laurent expansion as proposed in ref. || [7|. This 
expansion should be applied even in the case of amplitudes generated via the pinch 
technique. In ref.|| it was pointed out that the expansion was not just a mathemati- 
cal trick. The leading resonant term in the expansion represents the finite propagation 
and subsequent decay of a physical unstable particle. As such it constitutes a dis- 
tinguishable physical process and therefore must be exactly gauge- invariant. This 
insight was used to give an answer to the long-standing problem of how to calcu- 
late production cross-sections for unstable particles and was applied to the process 
e+e- -> Z°Z°. 

The presence of an unstable particle is indicated by a finite propagation length 
separating production and decay vertices. The contributions responsible for finite 
propagation are just those that are resonant in the invariant mass of the decay prod- 
ucts. Thus, for example, to calculate the cross-section e + e~ — > Z°Z° one calculates 
the matrix element for the process e + e~ — > Z°Z° — > (/i/i)(/2/2) and extracts the 
part resonant in the invariant masses, p\ and p\, of the (/i/i) and (/b/b) pairs respec- 
tively. This is accomplished by means of a Laurent expansion in each of the variables 
p\ and p\. In the notation of ref. the general matrix element for the process takes 
the form 

M(...,pi,p 2 , •••) = X^i(-,Pi,P2, ...)A l (t,u,pl,p 2 2 ) (1) 

i 

where the Zj are 'standard covariants', with external wavefunctions attached, that 
form a basis for the spinor and Lorentz tensor structure of the matrix element. They 
therefore transform as M. does. The Ai are Lorentz scalar functions of the indepen- 
dent Lorentz invariants of the problem and contain no kinematic singularities. The 
resonant part is extracted by taking the leading terms A, in a Laurent expansion first 
in p\ and then p\. 
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The cross-section formed from the resonant parts 

M(...,pi,pz, ...) = Y^li('..,Pi,P2, ...)Ai(t,u,p\,pl) (2) 

i 

or in this case doubly resonant parts, of the matrix element and then summed over all 
possible final states constitutes the cross-section a(e + e~ — ► Z°Z°). If the standard 
covariants U do form a basis then the scalars Ai will be gauge-invariant and hence so 
will the Ai. 

The external particle wave functions do not feature in the expansion process hence 
the kinematics of the problem are unchanged and momentum conservation is pre- 
served throughout. The final state integrations involve only stable on-shell particles 
and so things like complex scattering angles that have plagued certain other attempts 
at calculating such cross-sections do not arise. Problems near 'threshold' that have 
been noted elsewhere [0, do not appear. This is as expected since unstable par- 
ticles do not exhibit sharp thresholds as do stable ones. Branch cuts corresponding 
to their production lie off the real axis and such threshold problems must be calcu- 
lational artifacts. There is no ambiguity or flexibility at any stage of the procedure. 
Because matrix elements are always evaluated with external stable states they re- 
main gauge-invariant throughout. As guaranteed by Fredholm theory, there will be 
an exact factorization between the final state decay products of the unstable particle 
and the rest of the matrix element. Because of this exact factorization, the matrix 
element takes the form 



<t(s)= f dp\ f { S l] dp 2 2 a(s]pi,p 2 2 ) p(pl) p{jp\), (3) 

JO JO 

where cr(s;pi,p|) is the cross-section obtained from the resonant part of the matrix 
element alone and 

, ^ 1 p 2 (T z /M z ) 

p{p } ~ n-w^WTfm {po) {p ] 

is Breit-Wigner-like convolution function. 



In ref. [12] the process e + e~ — > u e e~W + — > u e e~ud is treated as an example of how 
to include width-effects in matrix elements involving unstable particles. Much is made 
of electromagnetic £/(l) e . m . gauge invariance and the possibility that the large mass 
ratio s/m 2 e present in this process might amplify uncancelled gauge- dependence to a 



disastrous level [13], [L4J]. In particular it is shown that gauge invariance guarantees 
that the cross-section behaves as ~ q~ 2 as q 2 — > and not ~ g~ 4 as it might otherwise 
do. However the question of how to simultaneously maintain SU{2) L and £7(l) e . m . 
gauge invariance in calculations that take account of W boson width effects is not 
considered. The approach advocated is to include a subset of higher-order corrections 



to restore U(l) c , m , gauge invariance as was done in ref.s[|l^, |lq| . This procedure is 
however inconvenient, difficult to apply consistently and adds considerably to the 
calculational labour involved. As noted in ref. the reparation scheme is arbitrary. 
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It may be commented that the use of next order diagrams to fix problems at a 
given order for processes involving unstable particles has been noted elsewhere fTF ]. In 
that case troublesome imaginary parts of counterterms were observed to be canceled 
by next order diagrams. 



2 Gauge Invariance for Massless Particles 

When massless gauge particles are present special problems arise. Suppose we were 
to write a Green's function G M with a photon leg in terms of some basis of standard 
covariants, If, and Lorentz scalars, A iy in a form analogous to eq.fll]) and 

G»( Pl ,p 2 ,...) = J2 1 ?(Pi,Pi>-)MpIpI-)- (4) 

i 

The If, of course, must transform in the same way that does. 

The gauge condition on Green's functions q^G^ 1 = means that the are not 
linearly independent since they satisfy the condition J2i(l ' h)Ai = where g M is the 
momentum of the photon leg. 

Gauge-invariant Green's functions can be obtained using the methods of Bardeen 
and Tung ||. Suppose that we have some exact Green's function expressed in terms 
of standard covariants and Lorentz scalars as in eq.(^). Since = the Green's 

function is invariant under of the action of the operator 

K = 9t-— (5) 
p-q 

where p is some conveniently chosen momentum. The spinor-Lorentz tensor structure 
of the Green's function is spanned by the set {I^lf}. That is 

G»( Pl ,p 2 , ...) = £ (I^( Pl ,p 2 , -..)) MpIpI •••)• (6) 

i 

For an approximate Green's function that does not satisfy the U (l) e . m . gauge-condition, 
the operator I£ serves as a projection operator onto {I„lf}- From the {I„lf} it is 
always possible to construct a new basis that is free of kinematic singularities that 
might have occurred when, for example, {p-q) =0. 

In general consider a Green's function G^ 1 calculated using some incomplete ex- 
pansion up to some order. It will consist of two parts 

G M = Go + Gf (7) 

where Gq is a consistent gauge-invariant contribution correct to given order of the 
calculation. Gf is a spurious higher-order gauge-dependent correction. Because Gq 
satisfies the gauge condition, q^G^ = 0, it is invariant under the action of the operator 
I". Thus at a given order we may make the replacement 

C -> I»G V = G% + I$G V X . 
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Only Gi is affected by the projection operator but since it is of higher order this 
is of no concern. The important point is that it now satisfies the gauge condition 
q^IffGi) = and is prevented from causing too strong a numerical blow up as 
q 2 - 0. 

3 The Process e + e~ —> u e e~W + — > u e e~ud 

Consider the process e + e~ — > v e e~W + — > v e e~ud. Some of the tree level diagrams 
contributing to this process are shown in Fig.l. In the kinematic regions where the 
process is allowed (but W + W~ production is forbidden) it is expected to be dominant 
because of the presence of a resonant W propagator. We will mainly concern ourselves 
with the region with q 2 small but the methods used apply in any kinematic region 
and will yield results that are exactly SU(2)l and £7(l) c . m . gauge-invariant order by 
order. 

Let q be the 4-momentum of the virtual photon and p + be that of the W + . The 
momenta of the incoming electron and positron are p\ and and k\ respectively and 
p 2 and k 2 are those of the outgoing electron and v e . Hence 

p+=p u +Pd, p- = h-k 2 , q = P\~P2- (8) 

M w is the W mass, e and g are the electromagnetic and SU(2) L weak coupling 
constants respectively. Qi is the electric charge of particle i. As usual jl — f (1 — 7s) 
is the left-handed helicity operator. 




Figure 1: Tree-level diagrams contributing to the resonant part of the process e + e — > 
v e e~W + — > v e e~ud. 

Using the same notation as in the introduction the exact matrix element takes 
the form 

M = Y,h(q,P+,...)A l (q 2 ,p 2 + ,...) (9) 

i 

This matrix element may be divided into a part that has a simple pole at q 2 = and 
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a part that has no pole 



These terms must be separately SU(2) L gauge-invariant because their differing pole 
structures. The first term receives contributions from diagrams, such as those of 
Fig.l, in which a photon is connected to the external electron current. The second 
term receives contributions from these as well many other topologies. Examples are 
the diagrams of Fig.l in which the photon is replaced by a Z° boson. 

We will be interested in the small q 2 region and will thus concentrate on the first 
term in eq. (|10|) and write 

m = Y,j;^r^p^-)R , l (o,P 2 + ,...) (ii) 

i 

where 

J7 = gr(0) u(P2)lMPi) / 12 x 

" [i-n 77 (o)]i' g 2 

The (^) M (g,p+, ••■) are standard covariants with attached external wavefunctions and 
Vy(q 2 ) is the 7e + e~ vertex form-factor. IT 77 (g 2 ) is the derivative of the photon self- 
energy with respect to q 2 . The factorization of the matrix element in this way is 
guaranteed by Fredholm theory. Since M. is £7(1) e .m. gauge invariant it is possible 
to find a set {(l' i ) fl (q,p + , ...)} that is free of kinematic singularities and that satisfy 
^)"(?>P + .-) = 0|. 

^From the part of the matrix element ([11]) we wish to extract the dominant reso- 
nant part of the matrix element that describes the production of a physical W + and 
its subsequent decay into the ud final state. As with the process e + e~ — > Z°Z° this 
is the leading term in a Laurent expansion of i?j(0,pi, ...) in p 2 ^, the invariant mass 
of the W + decay products. The result takes the form 

m = J2 j;j?(i"Y a (q,p + , -)R'm *w, •■•) (is) 



where 

jW = Vw(sw) u(PuhalLv(Pd) ^ 

[l-TV ww {s w )\^ p\ ~ s w 

Here U' ww (q 2 ) is the derivative of the W self energy with respect to q 2 and sw the 
solution of p\ — Myy — Uwwip+) — 0. Vw{p\) is W + ud vertex form factor. The 
result fljjD is exactly SU(2)l and U(l) e . m . gauge-invariant. 

Since the (l'-) fia (q,p + , ...) are U(l) e _ m , gauge-invariant they are invariant under the 
action of projection operators of the form given in eq.(|5|). 

In finite order calculations, because of the non-perturbative nature of the W + 
resonance, t/(l) e . m . gauge-invariance may not be satisfied. A gauge-invariant result 
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may be obtained by inserting a projection operator, Iff, as defined in eq.@. This 
has the the effect of discarding spurious gauge-dependent terms that must ultimately 
cancel in the final result thus avoiding the need to calculate higher-order diagrams as 
was done in ref.s[12, 15]. 

With these general considerations in mind we can now calculate the cross-section 
for the process a(e + e~ — > u e e~W + — > v e e~ud) in leading order and with the W width 
effects correctly taken into account. In words this is the production cross-section for 
e~z/ e and a physical W + that subsequently decays into a ud pair. 

As discussed above the contribution relevant as q 2 — > takes the form 

M = M^jp™. (15) 

in which 

2 

M ,KX = Y t Mt x - (16) 

i=i 

The Mf 01 correspond to the two diagrams of Fig.l, 

M^ a = ie^=Q w 2 1 V^ip^-p^-qMh^v^), (17) 
V 2 P- — J-VJ- W 

MT = te^QMh)Yj^^l a 7 L v(k 2 ). (18) 

(19) 

and 

V^ 2 ^( Pl ,p 2 ,p 3 ) = {pi-p 2 Y 3 g^ 2 + (P2-p 3 r i g^ 3 + (ps-pO'V 3 " 1 . (20) 

The projection operator I£ of eq.([|) may be inserted into eq.(^) in order to guarantee 
U(l) e .m. gauge invariance of the Green's function, M m . The matrix element then 
becomes 

M = M^I^JJJ^. (21) 
The matrix element squared and averaged over spins then takes the form 

\M\ 2 ) = \ £ M^Jlg apP {pl)M»WJj (22) 

e^e spins 



with 

2 ^2 



ud spins 



9 P+ „ n(^\n(J2 



48tt 2 \p\-s w \ 2 



g aP e(p° + )e(pl) (23) 



at leading order. 
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In forming the total production cross-section for a physical W + 's that subsequently 
decays into ud we must sum over the ud final states. Hence 

a(e + e- -> u e e~W + -> i/ e e~t«i) = N c J dp 2 a(p 2 )p(p 2 ). (24) 

The integral is over the full range allowed by the cut on the scattering angle, 9. N c 
is the number of QCD colours and 

2 If d 3 p 2 d 3 k 2 d 3 p + I 2 



a(j> ) 2s J (2vr) 3 2p° (2tt) 3 2£;° (2n) 3 2p° + V M ° 

x (2n) 4 5 4 ( Pl + h-p 2 -k 2 - p + ) (25) 

in which 



l-^o| 2 ) = i E M^llJlM^UJ^ (26) 



e 1 e spins 



and the integral is evaluated under the constraint that p\ = p 2 . Because Ad^I^ is 
^(l)e.m. gauge- invariant, (|.Mo| 2 ) ~ Q 2 only as q 2 — > 0. 
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